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SUMMARY 


The Limb Infrared Monitor of the Stratosphere (LIMS) is one of the experi- 
ments on the Nimbus-7 satellite. It is designed to scan the Earth's limb ver- 
tically and to measure spectral emission profiles of trace atmospheric gases 
that are believed to be important in processes controlling the stratospheric 
ozone distribution. The LIMS must have adequate spatial-frequency response for 
all the spectral channels to provide, through inversion of the measured limb 
radiance profiles, important information about the temperature, structure, and 
composition of the atmosphere. Experiment objectives are reviewed and several 
analyses and measurements are described which were performed to determine the 
adequacy of the system for satisfying these objectives. From the LIMS design- 
model data, the modulation transfer function (MTF) was calculated for the 
optical system, the detector field mask, the electronics, and the overall 
system for each channel. The signal output performance of the instrument was 
predicted from the system MTF data and model input radiance data for each chan- 
nel. The MTF measurements made on the flight sensor confirmed the analytical 
results. The predictions indicate that the instrument can satisfy the basic 
measurement objectives of the experiment. 


INTRODUCTION 

The experiment Limb Infrared Monitor of the Stratosphere (LIMS) , launched 
aboard the Nimbus-7 spacecraft in 1978, is designed to scan the Earth's limb 
vertically and to measure spectral emission profiles of trace atmospheric gases 
that are believed to be important in processes controlling the stratospheric 
ozone distribution. The limb radiance profiles of these gases (i.e., the 
spatial variation measured as the radiometer scans from the Earth's disk across 
the atmosphere) can be inverted (refs. 1 and 2 ) to provide important information 
about the temperature, structure, and composition of the atmosphere. The LIMS 
measures the radiance profiles of ozone (O 3 ) , water vapor (H 2 O) , nitrogen 
dioxide (NO 2 ) , nitric acid (HNO 3 ) , and carbon dioxide (CO 2 ) . The LIMS is one 
of several planned instruments (ref. 3) for surveying the atmosphere on a 
global scale and for obtaining an improved understanding of the physical and 
chemical relationships affecting its functions. The limited knowledge of the 
distributions of the gases within the atmosphere that LIMS will be measuring 
is due to natural variations such as latitude, season, and time, and to anthro- 
pogenic activities. For example, reference 4 provides a good survey of such 
variations for H 2 O. 

The specific experiment goals, system operation, and data-processing elec- 
tronics of LIMS have been described in reference 5. The present paper elab- 
orates on the spatial-frequency response characteristics of the LIMS and 
relates them to model radiance profiles of the type that the LIMS instrument 
is expected to observe. The limb-scanning technique used in the LIMS experi- 
ment, combined with the exponential falloff of atmospheric emission with alti- 
tude, provides an inherently high vertical-resolution potential (ref. 1 ). 



Design and performance trade-off typically includes field of view, signal to 
noise, target radiance, instrument collection aperture, and spectral range, as 
well as science requirements. 

After a brief review of the LIMS experiment objectives, the associated 
radiometric requirements and the analytical methods used to verify the optical- 
system design model are presented. Then analytical spatial-frequency response 
(modulation transfer function) data on the design model and the measured 
response data on the flight sensor are discussed. Finally, a brief assessment 
of the expected performance of the LIMS system is given, based on predicted 
inputs and on the system performance analysis. 

Throughout the LIMS system analysis, the authors relied upon the com- 
putational expertise of William L. Edmonds, Systems and Applied Sciences 
Corporation. 

Use of trade names or names of manufacturers in this report does not con- 
stitute an official endorsement of such products or manufacturers, either 
expressed or implied, by the National Aeronautics and Space Administration. 


LIMS EXPERIMENT OBJECTIVES 

Measuring the stratospheric ozone distribution is the primary objective 
of LIMS. A typical radiance profile for O 3 has been synthesized from studies 
reported in references 6 and 7 and from related studies supporting the develop- 
ment of LIMS and earlier limb-scanning experiments (refs. 8 and 9) . Typical 
radiance profiles have also been synthesized for H 2 O, HNO 3 , NO 2 , and CO 2 . The 
CO 2 radiance data, by virtue of the known uniform distribution of CO 2 , allow 
determination of the atmospheric vertical temperature profile needed for inver- 
sion computations. Carbon dioxide has been extensively modeled from temperature 
and pressure profiles representing temporal, geographical, and meteorological 
variations over the Earth's surface (refs. 6 and 7). These models account for 
the absorption of water vapor and ozone in the outer portions of the CO 2 spec- 
tral band. Representative radiance profile models for the gases that LIMS 
measures are shown in figure 1. For the purposes of this paper, the spatial- 
frequency contents of the radiance profiles are primarily of interest (i.e., 
their Fourier spectra, estimations of which are shown in fig. 2). The radiance 
profiles are synthesized from 64 data points which correspond to 1-km incre- 
ments of altitude; therefore, the maximum spatial frequency for which the 
spectra are determined is 0.5 cycle/km. It is seen that for all the profiles, 
the amplitude spectrum decays over the range defined by two or more orders of 
magnitude relative to the zero-frequency amplitude. The spectra were estimated 
by applying discrete Fourier transform techniques to the radiance profiles, 
after first multiplying the profiles by the Hann windowing function (ref. 10), 
followed by frequency averaging or smoothing (ref. 11 ). 

The LIMS experiment seeks to provide radiance measurements with a 12-bit 
(4096-level) amplitude resolution and with an error not greater than ±1 per- 
cent for constant radiance (dc) targets. In addition, the measurements should 
have sufficient resolution to follow atmospheric features with a period as 
short as 5 km in altitude (i.e., a spatial frequency of 0.2 cycle/km at the 
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atmospheric target for four narrow instantaneous fields of view (IFOV's) and 
0.1 cycle/km for two wide IFOV's). However, it is desirable that the instru- 
ment have the capability to respond to features of shorter period to determine 
whether they are present. 


LIMS RADIOMETRIC SYSTEM 

Figure 3 is a schematic of the LIMS optical system. The optical-system 
design was determined from a combination of resolution requirements, signal-to- 
noise constraints, and compatibility requirements set by the two-stage cryogenic 
cooler for the detectors (ref. 5). The system has three focal planes to satisfy 
the thermal and optical requirements. The first is located at the focus of the 
primary paraboloidal collector mirror, coincident with the chopper and calibra- 
tion source (not shown in fig. 3). The second focal plane coincides with a 
thermal mask that is important in the cryogenic cooling requirements. The final 
focal plane contains band- pass filters, the field-defining stop, and the detec- 
tors. The following optical- system studies that are described relate to the 
final focal plane. 

Figure 4 shows the shape of the field-stop mask and the corresponding 
detector array. The six detectors (four narrow and two wide) are located 
immediately behind the field stop. The detector dimensions are slightly larger 
than the field- stop apertures, so that all radiation within each IFOV is col- 
lected. The spectral band-pass filters are mounted on the front of the field 
stop. The atmospheric constituent corresponding to each detector is shown, as 
well as the wavelength value for each channel which was used in the optical- 
system analysis. Each wavelength in figure 4 is the midpoint of the LIMS 
channel spectral response. Because of the narrow spectral response of each 
LIMS channel, this midpoint corresponds essentially with the effective wave- 
length (centroid). Additional details of the LIMS optical system that are 
needed in the analyses are presented in appendix A. 

The LIMS will determine the vertical distributions of stratospheric O 3 , 

NO 2 , HNO 3 , H 2 O, and temperature with resolutions of about 1 km for temper- 
ature, O 3 , and HNO 3 , and of about 2 km for NO 2 and H 2 O. The CO 2 radiance pro- 
files for the two pass bands become data input to inversion algorithms to 
determine the vertical temperature profile of the atmosphere (refs. 1 and 2 ) . 

The inferred temperature profile is coupled with radiance data for the other 
gases to provide the distribution (density profile) of NO 2 , H 2 O, O 3 , and HNO 3 . 
The altitude ranges for measurements are from the troposphere to 75 km for 
temperature, to 64 km for O 3 , to 50 km for H 2 O, and to 40 km for HNO 3 and NO 2 . 


SPATIAL-FREQUENCY RESPONSE ANALYSIS 

The prediction and evaluation of electro-optical systems performance by 
optical transfer function (OTF) techniques are now recognized procedures 
(ref. 12). The OTF calculations (predictions) based on the design parameters 
were performed for the LIMS optical system with two different computer pro- 
grams for field angles corresponding to the geometric centers of the detec- 
tors. The two programs are POLYPAGOS (ref. 13) and ACCOS V (ref. 14). Before 
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comparative OTF calculations were done, agreement was established between the 
two programs for a paraxial ray-trace test case. For convenience, appendix B 
presents a brief tutorial on the analytical methods basic to these computer 
programs. 

The OTF calculations for the LIMS optical system were made for the direc- 
tion across the narrow dimension of the detector (i.e., along the vertical 
scan direction of LIMS) . The modulation transfer function (MTF) data for two 
typical channels (O 3 and H 2 O) of the LIMS that represent the two different 
IFOV's are shown in figures 5 and 6 . The agreement between the two sets of 
independent LIMS (optics) MTF calculations appears consisent with the current 
state of the computational art for complex optical systems. (See appendix B.) 
The phase transfer function (PTF) is shown for all the channels in figure 7. 
The MTF data are presented for the range of spatial frequencies of primary 
interest in the LIMS application. More complete data are included in appen- 
dix C. Two spatial-frequency scales, cycles/mrad and cycles/mm, are shown. 

The two scales are related by 


cycles 1 , cycles 

x - x 10 3 = ( 1 ) 

mrad f mm 


where f is the focal length (157 mm) of the optical system. Later these 
spatial-frequency scales will be related to the cycles/km (at target) , scale, 
in which performance requirements were specified. 

The MTF curve for the HNO 3 channel (appendix C) is nearly the same as 
that for O 3 . The MTF curves for the CO 2 channels (appendix C) are similar in 
shape to, but somewhat lower in magnitude than, the O 3 curves over most of the 
frequency range. The MTF curve for the NO 2 channel (appendix C) differs very 
little from that for the H 2 O channel over the frequencies of most significance. 

Appendix D presents methods used to define the best focal plane of the 
LIMS optical system. All the OTF data correspond to the best focus. 

The overall system response function for each channel is defined by the 
product of ( 1 ) the optical transfer function, ( 2 ) the response characteristic 
due to the geometric subtense of the IFOV stop, and (3) the response of the 
filter in the signal- processing electronics. Figure 8 shows the predicted LIMS 
total electro-optical system MTF for the O 3 channel. In this figure, the curve 
designated MTF(O) is the optical system MTF, taken from figure 5 (ACCOS V data). 
The MTF(S) curve represents the field-stop (detector-slit) contribution, given 
by the relationship (normalized) 


MTF (S) 


sin (7T$f) 
irBf 


(2) 
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where 3 is half the IFOV (mrad) and f is the spatial frequency 
( cycles/ mrad) . The curve designated MTF(F) represents the response of the 
four -pole Butterworth filter, which provides a cutoff (3-dB) frequency of 
4.4 Hz for the ozone channel, as well as for the CC >2 and HNO 3 channels. 

For convenience, the frequency scale is given in three different units. 
The temporal-frequency scale (Hz) is related to the angular spatial-frequency 
(cycles/mrad) scale by the relationship 


cycles/mrad cycles/mrad 

Hz = = (3) 

scan rate (mrad/sec) 4.36 mrad/sec 


The third frequency scale (cycles/km) is related to the cycles/mrad scale by 
the fact that 1 km at the target (Earth 1 s horizon) subtends 0.286 mrad at the 
radiometer for the nominal orbital altitude of 900 km 


cycles 


cycles x 0.286 mrad 


km (at target) mrad km (at target) 


(4) 


The MTF coordinate is shown over a range of three decades, 
are related by 


The two scales 


dB = 20 log (MTF) 


(5) 


Figure 9 presents corresponding predicted LIMS total-system data for the 
H 2 O channel. (See appendix C for similar results on the other four channels.) 
The detector- slit transfer function for the H 2 0 channel has a lower cutoff than 
that for the O 3 channel since the slit width for the H 2 0 channel is twice that 
for the O 3 channel. Also, the cutoff frequency of the electronics-filter func- 
tion for the H 2 O channel (NO 2 channel also (see appendix C) ) is at a value one- 
half (i.e., 2.2 Hz) that for the O 3 channel (CO 2 and HNO 3 channels also). The 
overall system response for the NHO 3 channel is practically identical with that 
for O 3 , and the system responses for the C0 2 channels differ insignificantly 
from the O 3 channel response. In addition, the NO 2 channel system response is 
practically identical with the H 2 O channel response. 

The dashed portions of the field-mask (slit) MTF curve (figs. 8 and 9) 
and of the system MTF curve represent the so-called "spurious resolution" 
regions of the system behavior beyond the first zero crossing of the field- 
mask transfer function. It is seen that the spurious resolution occurs at 
spatial frequencies greater than 0.57 cycle/km for the O 3 (CO 2 and HNO 3 ) 
channel and greater than 0.29 cycle/km for the H 2 O and NO 2 channels. The 
values of 0.57 and 0.29 cycle/km in the object (target) plane correspond to 
the values of 2 cycles/mrad and 1 cycle/mrad, respectively (i.e., the first 
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zero points in the slit transfer functions). The LIMS MTF (spurious) beyond 
these two spatial frequencies is down by at least three orders of magnitude 
relative to the zero-frequency MTF. 

The considerations in choice of the electronics-filter pass band are two- 
fold. On the one hand, the cutoff is chosen to pass all signal frequencies 
expected to be important. On the other hand, in the sampled data system used 
in LIMS (ref. 5) , the location of the 3-dB frequency, together with the steep- 
ness provided by four poles, mitigates any aliasing error, as illustrated for 
the O 3 channel in figure 10. The figure shows the overall MTF for the O 3 chan- 
nel folded about the sampling frequency of 41.7 Hz that is imposed by the LIMS 
Nimbus-7 data system. The logarithmic frequency scale distorts the mirror image. 
This figure indicates that aliasing errors are below the LIMS accuracy and 
resolution requirements since the overlap is nearly four orders of magnitude 
below the zero-frequency response. 

Figure 1 1 shows the phase response of the overall LIMS system for all the 
channels. These curves reflect mainly the phase response of the Butterworth 
filter since the PTF of the optical system itself (fig. 7) is nearly zero for 
most of the channels over most of the frequency range. 


MEASURED FREQUENCY RESPONSE 

The amplitude responses of the LIMS radiometric channels (flight sensor) 
were measured by using a series of input square waves of constant amplitude and 
differing frequencies. The test setup is sketched in figure 12. A series of 
thermal bar targets, consisting of alternating heated and cooled metal bars, 
was placed in the focal plane of the test collimator. The cooled portion of 
the target (back plate) was kept at liquid-nitrogen temperature to maintain the 
required constancy of input amplitude. The temperature of the warm portions 
(plate with slots) was controlled at 300 ± 1 K by means of a circulating bath. 
The measured data on the LIMS O 3 channel show some typical results. 

The response measured by means of the configuration shown in figure 12 is 
the square-wave response rather than the more conventional (sinusoidal input) 
MTF. The two responses can be compared by Fourier analysis. The response of 
the channel to the input square waves can be best understood by considering 
the following three frequency regions. 

(1) For square-wave input frequencies from 0 to about 0.01 cycle/mrad, the 
O 3 channel responds with little amplitude and phase distortion. Therefore, the 
output closely resembles the input square wave. Figure 8 shows that the third 
and fifth harmonics suffer hardly any decrease in amplitude. Figure 11 shows 
that the phase shifts for the third and fifth harmonics of the 0.01 cycle/mrad 
square wave with respect to the fundamental are about 13° and 25°, respectively. 

(2) For square waves in the range from 0.01 to 0.5 cycle/mrad, the output 
is distorted by increased phase shifts between the fundamental and harmonics. 

The phase shifts are caused mainly by the Butterworth filter. Phase shifts 
close to 180° between the fundamental and harmonics cause the resultant output 
wave to resemble a square wave with humps superimposed on it. The result is 
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that the maximum amplitude of the LIMS output for an input square wave in this 
frequency range is relatively larger than that for input waves in the lower 
frequency range, 

(3) For input frequencies greater than 0.5 cycle/mrad, the phase shift 
between the fundamental and harmonics decreases rapidly with increasing spatial 
frequency, but here the effects of steeply falling LIMS amplitude (sinusoidal) 
response suppress contributions of these harmonics almost completely. Conse- 
quently, for the higher frequencies, the output is a nearly sinusoidal response 
(i.e., a sinusoidal output is obtained for a square-wave input). 

The measured responses representing the three frequency regions previously 
discussed for the LIMS O 3 channel are shown in figure 13. The expected dis- 
tortions discussed are apparent in these measurements. The LIMS output response 
was characterized by determining the peak-to-peak amplitude of the output waves, 
which was then divided by the response for a zero-frequency target. These 
results are shown in table I for all the LIMS channels for spatial frequencies 
from 0 to 1.2 cycles/ mrad. 

The peak-to-peak measured square-wave response data for the ozone channel, 
corresponding to the data in table I, are shown in figure 1 4 as a series of 
circles. The solid-line curve, reproduced from figure 8 , is the MTF predicted 
for the O 3 channel from design data. The dashed curve represents the peak-to- 
peak square-wave amplitude response predicted from the ozone channel by the 
Fourier analysis considerations previously described. Peaking of the dashed 
curve in the midfrequency range occurs because of superposition and relative 
phase shifts of the square-wave harmonics. The measured response reflects a 
small increase in the midfrequency range, consistent with that predicted for 
square-wave inputs, and then a gradual decrease to sinusoidal-output behavior. 
However, the measured data (circles in fig. 14) are below the predicted square- 
wave response curve, consistent with the fact that system manufacturing and 
alignment errors would be expected to reduce the performance as compared with 
that of idealized predictions. The channel retains close to 30 percent of the 
zero-frequency response at the 1/(2 x IFOV) frequency (1 cycle/mrad) to reflect 
the proper trade-off with field-mask (slit) dimensions. 


PREDICTED FOURIER SPECTRUM OUTPUT AMPLITUDE 

To obtain the predicted output of LIMS for each channel, normalized to the 
zero-frequency output, the following relationship is used: 


Output spectrum = MTF (LIMS) x input spectrum 


( 6 ) 


This relationship applies because there is a spatially incoherent radiance 
input to the optical system (including field mask) which linearly maps object 
radiance into image irradiance at the detectors. The electronic system, in 
turn, linearly maps the irradiance into a voltage output, its transfer function 
being that of the Butterworth filter (figs. 8 and 9) . 
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Figure 15 presents the predicted normalized LIMS Fourier spectrum output 
amplitudes for all channels of the radiometer as determined from equation ( 6 ) 
(i.e.f from multiplying the calculated MTF by the predicted input spectrum 
amplitude for each channel (fig. 2)). These predictions indicate that, for 
the H 2 O and NO 2 channels, the output spectrum magnitude (relative to the zero- 
frequency output) at 0.286 cycle/km, where the resolution becomes spurious, 
would be less than 10"^. For the O 3 and HNO 3 channels, the output spectrum 
magnitude at 0.572 cycle/km, where the resolution becomes spurious, would be 
about 2 x 10 ~ 5 and 10 -5 , respectively. The CO 2 channel output at 0.572 cycle/km 
appears to be less than 10 - ® because its spectral falloff is faster than that 
for the other narrow IPOV channels. These predicted outputs are based on some 
extrapolation beyond the spectral range of the input- spectrum calculations, as 
indicated by the dashed parts of the curves in figure 15. The implications are 
that any spurious resolution would be indistinguishable from noise unless the 
system could accurately quantize more than 50 000 radiance levels. 


As stated earlier, the LIMS experiment seeks to provide radiance- 
measurement capability, with an amplitude resolution of 12 bits, at spatial 
frequencies at the atmospheric target up to 0.2 cycle/km for the narrow IFOV' s 
and up to 0.1 cycle/km for the wide IFOV's. Nominal variations of the instru- 
ment temperature in orbit cause changes in the self-radiation of the optical 
system, to which the infrared detectors are sensitive. To prevent possible 
saturation of the channels due to variations in the optical-system radiation 
component, the full-scale radiance excursion (midrange) of each channel is set 
to 80 percent of the full 12-bit (4096-level) data-system range or to 3276 quan- 
tization levels (counts). Therefore, the operational requirement on amplitude 
resolution becomes one part in 3276, or 3.1 x 10 -4 , up to 0.10 cycle/km for the 
NO 2 and H 2 O channels, and up to 0.20 cycle/km for the O 3 , CO 2 , and HNO 3 chan- 
nels. As seen in figure 15, the predicted LIMS output magnitudes at these 
spatial frequencies, normalized to zero-frequency output, are 5.4 x 1 0~2 f or 
NO 2 , 4.3 x 10~2 for H 2 O, 2.6 x 1 0 - ^ f or 03 , 4.1 x 1 0 -4 for CO 2 , and 3.1 x 1 O - ^ 
for HNO 3 . Therefore, the predicted outputs for all the channels are greater 
than the required amplitude resolution at the highest spatial frequencies of 
interest by factors ranging from slightly more than 1 for the CO 2 channel to 
174 for the NO 2 channel. These predicted outputs indicate the LIMS performance 
is adequate, provided the noise levels for all the channels are also less than 
the required amplitude resolution. Data that have become available since the 
October 1978 launch of Nimbus-7 (personal communication with John C. Gille, 
National Center for Atmospheric Research, and James M. Russell III, Langley 
Research Center, who were principal co- investigators for the LIMS experiment) 
show the typical noise levels relative to the maximum channel range are 
3.8 x 1 0 -4 for NO 2 ; 4.2 x 1 0 -4 for H 2 O; 4.6 x 1 0 -5 for O 3 ; 8.9 x 1 0 “ 5 for 
C 02 (wide); 3.6 x 1 0 -4 for CO 2 (narrow); and 7.0 x 10 -5 for HNO 3 . The noise 
levels represent the total spatial-frequency bandwidth of each channel. 

The noise levels appropriate to the 0.10 cycle/km (NO 2 and H 2 O) and the 
0.20 cycle/km (O 3 , CO 2 , and HNO 3 ) frequencies would be much further reduced 
by the effects of the electronics MTF on the noise spectrum. 
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CONCLUDING REMARKS 


Experiment objectives of the Limb Infrared Monitor of the Stratosphere 
have been reviewed and related to its predicted and measured spatial-frequency 
response. Consideration of the predicted output signal spectrum and the 
flight-observed noise levels together indicates very favorable signal- to- noise 
characteristics up to the highest significant spatial frequencies. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
June 8, 1979 
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APPENDIX A 


LIMS OPTICAL SYSTEM 

Figure AT repeats the LIMS optical- system schematic shown in figure 3 but 
breaks it into three parts to show all the important construction parameters 
needed for the optical transfer function (OTF) calculations. Figure AT (a) 
shows the front-end elements, including the first lens, and figure A1 (b) shows 
the remainder of the system, with the first lens also included. Figure A1 (c) 
shows details of the last lens in the system and of the focal-plane assembly 
(i.e., the spectral filters, field mask, and detectors). The field mask that 
defines the instantaneous field of view (IFOV) for each channel was shown in 
figure 4. The view of the field mask shown is along the optical axis Z 
(i.e., perpendicular to the view in fig. A1 (c) ) . The field-mask configura- 
tion is symmetric about the X- and Y-axes. 

Table AI shows the refractive indexes at the appropriate temperatures 
(see fig. 3) and wavelengths for the CdTe lenses and window and for the Ge sub- 
trates of the spectral band-pass filters. These index data, along with the 
temperature data in figure 3 and the information in figures 4 and Al , provide 
all the necessary inputs for calculating the OTF of all channels of LIMS. 


TABLE AI.- INDEX OF REFRACTION 


Wavelength, 

Vim 

Index of refraction for - 

CdTe (Irtran-6) 1 at 

temperature, K, of - 

297.5 

152 

65 

Ge at 
65 K 

6.21 (N0 2 ) 

2.68499 

2.671 1 5 

2.66288 

3.91 743 

6.83 (H 2 0) 

2.68332 

2.669533 

2.661 294 

3.915565 

9.64 (0 3 ) 

2.67621 

2.66236 

2.65409 

3.91 131 

H.32 (HN0 3 ) 

2.67139 

2.65755 

2.64927 

3.91028 

15.20 (C0 2 wide) 

2.65751 

2.64390 

2.63575 

3.90934 

15.23 (C0 2 narrow) 

2.65751 

2.64390 

2.63575 

3.90934 


^Irtran-6: Registered trademark of Eastman Kodak Company. 
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View 

before 

Z-axis 

folding 

mirror 




(a) Input end. 

Figure A1 LIMS optical- system schematic. All dimensions are given in centimeters. 
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- Vertex 
radius of 
curvature: 
12 . 1923 - 



(b) Output end. 


Figure A1 .- Continued. 
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Radius of 
curvature: 
.9157- 


Radius of 
curvature: 
.5977 


— .2537 


.4191 



Detector 
plane - 



22367 (H 2 0) 


21986 (0 3 ) 


22225 (HNO a ) 


.23569 (C0 2 (Wide))- 
23569 (C0 2 (Narrow)) 


.3429 


.22200 (N0 2 



Filters 
(Ge) 

Filter 
thickness 


Field -of -view 
mask (best-focus 
plane) 


Filter 

band 

Filter 

thickness 

H 2 ° 

0.07453 

°3 

.07833 

hno 3 

.07595 

C0 2 (Wide) 

.06251 

C0 2 (Narrow) 

.06251 

no 2 

.07620 


(c) Focal-plane assembly 
Figure AT Concluded. 
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ANALYTICAL METHODS FOR OPTICAL TRANSFER FUNCTION CALCULATIONS 

Although the calculation of the optical transfer function (OTF) from 
optical design parameters is, in principle, a straightforward matter, many 
approximations and simplifying assumptions are required, in practice, to make 
the numerical calculations manageable. Different analysts and designers make 
different compromises in their computer programs and generally get somewhat 
different results. Only recently have these differences been well enough 
understood to provide good agreement between results obtained by different 
evaluators and to assure good reproducibility among results of an individual 
evaluator (refs. 12, 15, and 16). 


POLYPAGOS Computer Program 

POLYPAGOS (currently named the Spectral Lens Analysis Program (SLAP) and 
developed by the Pagos Corporation) is a general-purpose program for analysis 
of complex optical systems. Specific mathematical details can be found in 
references 13, 17, and 18, and background treatments can be found in such text- 
books as references 19 and 20. With POLYPAGOS and similar programs, the first 
step in the analysis of a system is the computation of the pupil function by 
suitable ray-tracing procedures. The pupil function describes the combined 
effects of diffraction (system limiting aperture) and aberrations on the other- 
wise ideal wavefront propagated through the system. 

Fourier transformation of the pupil function yields the amplitude spread 
function (ASF) , and the point spread function (PSF) is obtained from 
PSF = | ASF| 2 . Physically, the PSF is the diffraction pattern or optical power 
spectrum in the image plane for a point-source object (input) , as would be 
measured with a square-law detector such as photographic film or a photo- 
detector. The OTF which results from a Fourier transformation of the PSF is 
a two-dimensional complex function that may be written 


OTF = R(f x ,f y ) + il(f x ,f y ) (Bl) 

where R and I represent the real and imaginary parts and f x and f y are 
spatial-frequency variables for the x- and y-directions. The OTF must be 
specified for a given field position, for a given direction, and for a given 
wavelength unless a polychromatic quantity is determined by averaging. 

The phase transfer function (PTF) is given by 


4 > 


tan -1 


I(f x ,f y ) 

R(f x »fy) 


(B2) 
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and the modulation transfer function (MTF) by 
MTF = [R2(f x ,f y ) + l2(f x ,f y ^| 1/2 


(B3) 


The meanings of MTF and PTF are best understood by considering that a given 
distribution of radiance in the object plane of an optical system can be con- 
sidered as arising from the superposition of sine-wave grating distributions of 
radiance of different spatial frequencies, phases, orientations, and amplitudes. 
The distribution of irradiance in the image can be built up from sinusoids in 
the same way. The MTF is the variation with spatial frequency of the ratio of 
the contrast in the sinusoidal image that is produced to the contrast of a sinu- 
soidal test object. 

The PTF is a measure of the lateral displacement of the component sinusoids 
from the origin of the image coordinate system; it is measured in radians or 
other angular units, where 2^ radians are equivalent to one period of the image 
of the sinusoidal test object. The PTF is the least understood aspect of the 
OTF. Reference 15 (pages 39-43) shows that the meaning of the PTF is always 
ambiguous unless the PSF is evenly symmetrical. Therefore, the practical impor- 
tance of the PTF in image evaluation is still to be fully resolved. A linear 
phase shift just translates the diffraction image without detrimental effect 
on the image quality. Programs such as POLYPAGOS provide for automatic removal 
of linear phase shifts. Nonlinear phase shifts, however, introduce harmonic 
distortion that results in serious deterioration of image quality. 


ACCOS V Computer Program 

The ACCOS V program is a comprehensive program including optimization rou- 
tines for use in designing and evaluating a broad range of optical systems. 

The ACCOS V program derives the OTF by autocorrelation of the pupil function, 
a process mathematically equivalent to the Fourier transformation operations 
of the POLYPAGOS program. However, if the full two-dimensional OTF is not 
required, the autocorrelation process is more computationally efficient than 
the Fourier transform process. 
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APPENDIX C 


SUPPLEMENTARY MODULATION TRANSFER FUNCTION DATA 
ON LIMS RADIOMETRIC SYSTEM 

In the main text of this report, modulation transfer function (MTF) data 
were presented on the O 3 and H 2 O LIMS channels over the spatial-frequency range 
of primary significance in the LIMS application (figs. 5 and 6 ) , and the close 
similarities of data on the other channels to these two data sets were dis- 
cussed. In this appendix, additional MTF data for all six channels are pre- 
sented for spatial frequencies up to about the cutoff spatial frequency of the 
optical system for the CO 2 wavelength. The data might be of interest (perhaps 
aside from the LIMS application) in the frequency response of a LIMS-type 
radiometric system. These results are contained in figures Cl to C 6 . 

In addition to calculations of MTF at field angles corresponding to centers 
of the detectors (represented by the data points in figs. Cl to C 6 ) , MTF values 
were calculated for field angles corresponding with the detector extremities 
(see circular dots on detectors in fig. 4) ; the extremity values of MTF are 
represented in figures Cl to C 6 by the ends of the vertical bars. The spread 
(y-direction) in MTF values for different field angles along the length 
(x-direction) of each detector is due to the variation in aberrations with 
increased or decreased field angle. 

The differences between the POLYPAQOS and ACCOS V MTF results may be 
accounted for through the considerations of the first paragraph of appendix B. 
However, the differences are not very significant relative to the final system 
MTF results (i.e., the MTF (LIMS) curves in figs. 8 and 9 and figs. C7 to CIO). 
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Figure C3.- MTF(optics) of LIMS O 3 channel. 
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Figure C4.- MTF (optics) of LIMS HNO3 channel. 
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Figure C5.- MTF (optics) of LIMS CO2 (wide) channel. 
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APENDIX D 


OPTIMIZATION OF FOCAL PLANE FOR OPTICAL 
TRANSFER FUNCTION CALCULATIONS 

All the optical transfer function (OTF) calculations on the LIMS optical 
system that are presented in figures 5 and 6 and figures Cl to C6 were obtained 
for the best focus. The best focus is defined as the focal plane for which the 
modulation transfer function (MTF) data for all the detectors are most nearly 
maximum for a selected spatial frequency. The spatial frequency chosen was that 
equal to the reciprocal of twice the detector mask width (y dimension) . For the 
NO 2 and h 2° detectors, the reference spatial frequency is 0.5 cycle/mrad, and 
for the other four detectors it is 1 cycle/mrad. For a given detector, this 
spatial frequency is one-half the first cutoff in the field-mask (slit) trans- 
fer function. (See figs. 8 and 9.) 

Figure D1 shows the MTF, calculated for the reference spatial frequencies 
previously cited, as a function of focal position (z-direction) for all the 
detectors, as determined with the POLYPAGOS program and verified with the 
ACCOS V program. It can be seen that the best-focus choice (zero z position, 
abscissa of fig. D1 ) is a slight compromise among the six detectors. However, 
the 20- to 30-Mm variation in best focus among the six detectors represents 
only about 2 units of variation in MTF (i.e., 2 percent of the zero-frequency 
value) . 

Figure D2 presents a further verification of the best focus choice, which 
is determined by using a POLYPAGOS subroutine to calculate the percentage of the 
radiant power from a point-source object, in an appropriate field position, 
that is available at each detector (ordinate) . The zero z point (abscissa) 
in this figure is chosen identical to that in figure D1 . The peaks on all these 
curves lie within 20 jim of the chosen best focus, and the peak MTF value for 
any curve is not different from the value at z = 0 by more than 1 percent of 
the zero-frequency MTF value. Note that the ordinate scale in figure D2 is 
expanded relative to that in figure D1 . 
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Figure 3.- LIMS optical schematic. 







4.1 mrad (.6425 mm) 




Figure 4.- Detector (f ield-of-view) mask and detectors. The black dots 
represent field points for OTF calculations and X represents central 
wavelengths. 
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Figure 5.- MTF (optics) of LIMS O 3 channel. 






PTF 




39 


I 



0 


H 

S 


I 


0 


.01 


L .001 


.01 


□ MTF(0) 

A MTF(S) 

O MTF(F) 

X MTF(LIMS) 


i i i m l 

.i 




Mini i i 1 M i ; il I II 1 1 ijjlL 

.1 1 10 


10 


Cycles 

mrad 

Hertz 


i ii 1 1 ill i i i i lull! i j _i l.i nil i 

.01 .1 1 

Frequency 


Cycles/km 
(at target) 


Figure 8.- Frequency response of LIMS radiometer for O3 channel. 
(Desired pass band is 0.0 to 0.2 cycle/km.) Dashed curve 
indicates spurious resolution. 
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Figure 9.- Frequency response of LIMS radiometer for H 2 O channel. 
(Desired pass band is 0.0 to 0.1 cycle/km.) Dashed curve 
indicates spurious resolution. 
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Figure 11.- Phase transfer function of LIMS radiometer. 
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Figure 12.- MTF test configuration. 
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(a) Target segment 1 . 

Figure 13.- LIMS output of MTF test 
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Figure 14.- Frequency response of LIMS radiometer for O 3 channel. 
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Figure 15.- Prediction of LIMS spectral output. Dashed line 
indicates extrapolation of data. 
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